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1. INTRODUCTION 

Information about the internal structure of the nucleon is encoded in its structure functions. 
While they cannot be computed directly on the lattice, the operator product expansion (OPE) pro- 
vides a connection between their moments and nucleon matrix elements of local operators. For 
instance, for the unpolarised structure function F\ , the OPE reads 

2 jxx"-ifi(^,e2) = £4f„vl^)+o(i/e'), (1.1) 

if) 

where / denotes the quark flavour, Ep^ ^ is the (perturbative) Wilson coefficient and the matrix 

If) 

element v,, is defined by 

(A^(p)|o|j')'-^"^ -traces|A^(^)) = 2v[!\p^' . . . p^" - txnct?.) , (1.2) 

with the operator 

OjJ-"''" = i/iyy^'i D^^ ... D^" xj/f. (1.3) 

Similar relations hold for the other structure functions, see e.g. M] for details. Both the matrix 

If) (f) 
element v,Y and the Wilson coefficient Ep^ ^ depend upon the choice of a renormalisation scheme 

and scale; only in their product, these dependencies cancel. 



2. LATTICE SIMULATION 

We use the overlap operator given by 



= P(l + ^ + (l- ^)75Sgn(//w(-p))) , 



(2.1) 



where Hy^{—p) = 75 (^^w — p), Dw being the Wilson Dirac operator. We approximate the sign 
function appearing in ( |2. 1[ ) by minmax polynomials [^. For the gauge part we chose the LUscher- 
Weisz action ^ 



S[U] 



CO £ iReTr [1 - U^u^] + ci £ iReTr [1 - U,,,,] + C2 £ iReTr [1 - U^.^] ] , (2.2) 

^ plaq rect 



pai' 



with coefficients ci, C2 (cq + 8ci + 8c2 = 1) taken from tadpole improved perturbation theory 
We ran our computations at two lattice spacings, see table |l|. The scale has been set from the pion 
decay constant, details are given in [@]. The parameter p in (0|) was set to 1.4. 



V 


i8 


a(fm) 


confs. 


16^32 
24^48 


8.00 
8.45 


0.153(3)fm 
0.105(2)fm 


300 
200 



Table 1: Parameters of the gauge configurations used. 
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In order to remove 0{a) errors from the three-point functions, we employ the method of 
which amounts to replacing propagators D^'*F by j^-nrD^^^ — ni^^m Jacobi smeared point 
sources [|I|] with parameters Ng = 50 and = 0.21 have been used in order to obtain a good overlap 
with the ground state. 

The computation of matrix elements follows the procedure outlined in ||l|] : We form the ratio 

^ ^ {N{t,ink)0{T)N{t,^ource)) 
{N{t,ink)N {Source)) ' 

from which the matrix element can be extracted in the region t source < f < tsi„u. We always set 
t source = and tsink = 9 {tsink = 13) on the j8 = 8.0 (j8 = 8.45) configurations (in lattice units), 
which corresponds to a distance between source and sink of 1 .4 fm. 

The matrix elements we are considering are listed in table along with the operators used for 
their determination, where we use the operators (L3) and 

We compute only flavour non-singlet matrix elements, since in this case there is no contribution 
from disconnected diagrams. 



3. NON-PERTURBATIVE RENORMALISATION 



The operators appearing inside the three-point functions have to be renormalised. For gA = 
Au — Ad, the operator to be used is the axial current , the renormalisation of which is particularly 
simple because it does not depend on a renormalisation scheme or scale. It can be obtained from a 
Ward identity [^] as 

2m, (P(OP(0)) 



lim lim 



(3.1) 



The renormalisation constants of the other operators under consideration are logarithmically 
divergent. We have computed them in the RI' — MOM-scheme [||]. In this scheme, the renormali- 
sation condition is formulated in terms of quark Greens functions in Landau gauge with an operator 
insertion at zero momentum transfer: 



Coip) 



-ip{x-y) 



{Yix)Oizmy)). 



(3.2) 



Matrix Element 


Operator 


gA 




gT 




V2 




a\ 


i(or+of) 



Table 2: Operators used in the three-point functions 
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(GeV ) 



ti'(GeV-) 



Figure 1: Left: The renormalisation constant Z^^ ^'^^{n). Right: The same, but with the renormahsation 



group running removed according to (3.6) 



From this quantity, the amputated vertex function Tq is formed: 

To{p)=S-\p)Co{p)S-\p), (3.3) 

with the quark propagator 



s{p) = -l^e-^''^^-y\wm{y))- (3.4) 



The renormalisation condition at scale }JL is 



Z^(At)Zo(At)no(ro(p)) 1^2=^2 = 1, (3.5) 



with the projector no(r) = jj}^ \^ o^Bom^P)^)- "^^'^ wavefunction renormalisation constant Z,^ 
has been determined from the relation Z^,ZaYIa (F^) = 1. 

In order to convert the results to the MS-scheme, we first determine the renormalisation group 



invariant renormalisation constant Z^^: 



Zg^i = (zf -™''^^i(At))"'z^'-™(M), (3.6) 



and then convert to the MS-scheme at scale jx' (we always use fx' = 2GeV): 
ZMS(^') = Z^^''*°^(At)Zg°i with the conversion functions 



/ 

Jo 



(3.7) 



The coefficients of the j8 and / functions are taken from ||9|,[T0|]. 

In fig. [I] we plot zf^^'^^^ and Zf^^ for jS = 8.45. From the latter plot, we read off Zf^^ = 
2.6 from the plateau region SGeV^ < p^ < 15GeV2. Using zjff ''*°\2GeV) = 0.737, we obtain 
Z^^ = 1.92. The renormalisation constants for all operators we need are shown in table |[ A 
comparison with results obtained in one-loop tadpole-improved lattice perturbation theory [ pi] , [T2| ] 
shows large discrepancies, especially for the operators with one derivative. 
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O P=8.0 
□ R=8.45 



m"(GeV)" 



H-0 
□ |3=8.45 



m (GeV)" 



Figure 2: The nucleoli's axial charge (left plot) and tensor charge (right plot) as a function of the squared 
pion mass 



4. RESULTS 

Our results for the axial charge gA are displayed in fig. ^ (left). A linear extrapolation to the 
chiral limit yields gA = 1.37(5) at j3 = 8.0 and gA = 1.14(5) at j8 = 8.45. The tensor charge gj = 
5u — 5d is plotted in fig. ^ (right); in this case, Unear extrapolations to = yield gr = 1.35(4) 
at 15= 8.0 and gr = 1.18(5) at jS = 8.45. 

The results for the matrix elements Vj^^ = (x)"^^ and a'p^ = 2{x)^"^^'' are shown in fig. ||, 
together with the phenomenological values. In both cases, there is almost no dependence on the 
quark mass visible. In the range rrij^ > 400MeV, our results agree with previous results obtained 
from improved Wilson fermions 



5. CONCLUSIONS 

We have determined the flavour non-singlet nucleon matrix elements gA, gr, ^2 and a\ from 
quenched overlap fermions. The renormalisation has been done nonperturbatively. Comparing the 
results at the two values for the lattice spacing we have, we find significant discretisation effects, 
in contrast with the situation for hadron masses 

While our results are in good agreement with earlier determinations, there remains a rather 



Operator 


i8 = 


= 8.0 


i8 = 


8.45 


pert. 


nonp. 


pert. 


nonp. 




1.36 


1.59 


1.30 


1.42 




1.36 


1.73 


1.33 


1.54 




1.33 


2.11 


1.39 


1.92 




1.34 


2.21 


1.40 


1.98 



Table 3: Comparison between the renormalisation constants Z^^(2GeV) obtained non-perturbatively and 
in lattice perturbation theory. 
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O P=8.0 



|3=8.0 
□ P=8.45 



m^{GeV) 



m7GeV) 



Figure 3: The nucleon matrix elements a\ ^ (left plot) and (right plot) as functions of the squared pion 
mass (MS, 2GeV). 



large discrepancy to the phenomenological values for V2 and ai, even at the lowest quark masses 
we can reach at present. 
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